Minimum energy states of the plasma pinch in standard and Hall 

magnetohydrodynamics 
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Axisymmetric relaxed states of a cylindrical plasma column are found analytically in both stan- 
dard and Hall magnetohydrodynamics (MHD) by complete minimization of energy with constraints 
imposed by invariants inherent in corresponding models. It is shown that the relaxed state in 
Hall MHD is the force-free magnetic field with uniform axial flow and/or rigid azimuthal rotation. 
The relaxed states in standard MHD are more complex due to the coupling between velocity and 
magnetic field. Application of these states for reversed-field pinches (RFP) is discussed. 



Minimum energy states appear in many magnetized 
plasmas as a result of relaxation - the spontaneous ten- 
dency to evolve toward preferred configurations with or- 
dered structure. Theoretical prediction of such states is 
a long-standing problem for both laboratory and astro- 
physical applications. The concept of relaxation was pro- 
posed by Taylor [1], who conjectured that during turbu- 
lent dynamics a slightly resistive magnctohydrodynamic 
(MHD) system tends to minimize its magnetic energy 
while conserving the total magnetic helicity. The under- 
lying basis of this approach is the principle of selective 
decay of invariants [1-3], i.e., one or more ideal invariants 
of the system (conserved in the absence of dissipation) 
are less susceptible to dissipation than energy and thus 
can be considered as constants during the relaxation pro- 
cess. Mathematically the relaxation theory is formulated 
as a variational procedure for obtaining a relaxed state 
by minimizing the energy subject to constraints. 

The Taylor theory is successful in explaining magnetic 
structures in laboratory plasmas such as the reversed- 
field pinch (RFP), multipinch and spheromak [2, 3] but 
it does not include flows that are ubiquitous in exper- 
imentally observed relaxed states. The origin of these 
flows is not well understood; laboratory plasmas rotate 
in the toroidal and poloidal directions even in the ab- 
sence of externally applied torques (intrinsic rotation). 
Further, the experimental parameters are such that the 
single-fluid MHD model may not be strictly valid, and 
the inclusion of the effects of separate ion and electron 
fluids in the model may be required. 

The present work is motivated by the recent progress in 
plasma velocity measurements in the Madison Symmet- 
ric Torus RFP experiment, which show an abrupt change 
of the global flows during the relaxation events. De- 
tailed temporal and spatial measurements of flow dynam- 
ics indicate significant radial angular momentum trans- 
port and flattening of the radial flow profiles [4]. 

The goal of the present paper is to determine the min- 
imum energy (relaxed) states for a cylindrical RFP, to 
analyze the possibility of plasma flows in such states in 
both standard (single-fluid) and Hall MHD (a two-fluid 



model with massless electrons and uniform density), and 
to elucidate their global properties. Since the RFP has 
"bad" magnetic curvature everywhere, the geometry of 
periodic cylinder is a good approximation for RFP theory 
and simulations [2, 3] . We employ a variational procedure 
that includes all ideal invariants inherent in correspond- 
ing models. While the experiments are open systems that 
interact with the external environment through applied 
voltages, here we consider only closed systems. This is 
consistent with Taylor's approach [1, 2], which has been 
reliable for predicting the general properties of relaxed 
magnetic fields without flow. The fields and flows pre- 
dicted by the present theory may be relevant to the flows 
that are observed after the crash phase of sawtooth cycle 
in the RFP. Of course, we cannot comment on the ori- 
gin of these flows, only on their relaxed properties. For 
simplicity, only axisymmetric states are considered. 

In the framework of incompressible MHD similar stud- 
ies are reported in Refs. [5-7]. In Ref. [5] (results are 
corrected in Ref. [6]) the cross helicity invariant is in- 
cluded in the analysis and a relaxed state with veloc- 
ity parallel to force-free magnetic field is obtained. We 
emphasize here that the cross helicity in incompressible 
MHD is an ideal invariant of the system and is also a 
rugged invariant in the presence of dissipation, which is 
confirmed by numerical simulations [8] . In Ref. [7] in ad- 
dition to the cross helicity the total angular momentum is 
taken into account as a conserved quantity in toroidal ge- 
ometry. This yields two separate types of relaxed MHD 
states with force-free magnetic field and either parallel 
flow or rigid rotation. The novel result of our paper is 
that the inclusion of the cross helicity along with mo- 
menta constraints in cylindrical geometry generally leads 
to MHD relaxed state with non force-free magnetic field 
and mixed types of flows. 

The relaxation problem in the framework of incom- 
pressible Hall MHD is considered in Refs. [9-13]. In 
Refs. [9, 10] the doublc-Beltrami equation for the relaxed 
states in Hall MHD is revealed and several solutions are 
presented for cylindrical geometry. These solutions do 
not correspond to true minimum energy states for fixed 
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electron and ion helicities, H e and Hi, since the mini- 
mization procedure is not completed. This is because 
the unknown Lagrange multipliers used in the variational 
principle are not specified in terms of the initial values of 
the invariants. A more complete analysis is reported in 
Ref. [11] where the energy of relaxed states is found in 
toroidal systems as a function of H e and Hi. Although 
the general solution of the double-Beltrami equation has 
two eigenfunctions, Ref. [11] uses only one of them. This 
precludes the possibility of two different spatial scales in 
the relaxed state (as in Refs. [10, 12, 13]). In Ref. [12] 
relaxed states are obtained as a linear combination of two 
orthogonal Beltrami eigenfunctions with eigenvalues Ai 
and A2, respectively, and the energy is expressed as a 
function of Hi, H e , Ai and A2. The next step is to find a 
pair of eigenvalues (Ai, A2) that minimizes the energy; 
however, this step is missing in Ref. [12]. 

The full energy minimization of the Hall MHD sys- 
tem with fixed electron and ion helicities is completed in 
Ref. [13], with the result that the relaxed state in Hall 
MHD is always a force-free magnetic field with no plasma 
flows. The authors of Ref. [13] question the conservation 
of the ion helicity Hi, arguing that it is not a rugged 
invariant during relaxation and, therefore, it should not 
be included into energy minimization procedure. The 
fact that it is not conserved in Hall MHD relaxation is 
confirmed by numerical simulations [13-15]. We adopt 
these conclusions by excluding the ion helicity from the 
invariants of the Hall MHD relaxation. Instead, we in- 
clude other velocity related invariants that follow from 
the geometrical symmetry of the system: total axial and 
angular momenta. This allows us to obtain a relaxed Hall 
state with plasma flows. 

We consider the problem of finding the minimum en- 
ergy (relaxed) states of cylindrical plasma pinch of length 
L. We assume that plasma is surrounded by a perfectly 
conducting shell (flux conscrver) of radius a, and plasma 
density p is uniform in space and constant in time. Un- 
der these assumptions the plasma is described by equa- 
tions of ideal incompressible Hall MHD, which in non- 
dimensional form are [9] 

dv 

— + (v-V)v + Vp=(Vxb)xb, V-v = 0, (1) 

OT 

— = Vx(vxb - 6j(V x b) x b), V ■ b = 0. (2) 
Here, all physical quantities are normalized: 



with $0 being a total axial magnetic flux. System (l)-(2) 
also contains the non-dimensional ion skin depth (or Hall 
parameter) Si, which is defined as 



s = <k = c / rr?j 
1 a ay Airpe 2 Z 2 ' 



where m, and eZ are ion mass and charge, c is speed of 
light. In the case Si = 0, single-fluid MHD is recovered. 

We adopt cylindrical coordinate system {r, tp, z}. In 
order to solve Eqs. (1), (2) uniquely, we have to specify 
boundary conditions for velocity and magnetic field. At 
a perfectly conducting, impermeable surface the normal 
components of the velocity and the magnetic field vanish, 



0, b n = 0, 



(4) 



and for closed systems the tangential component of elec- 
tric field is zero, which is equivalent to 



htin x b t ) = 



Sijn = 0, 



(5) 



where j = V x b is normalized current density. Note that 
in single-fluid MHD (Si = 0) this condition is satisfied 
automatically. 

Following the fundamental idea of the Taylor theory, 
we introduce a number of global invariants and examine 
them for cylindrical pinch geometry assuming boundary 
conditions (4), (5). The total energy of system (1), (2), 



E 



(v 2 +b 2 )d 3 r, 



(6) 



is an ideal invariant. In a dissipative isolated plasma, 
total energy can only decrease in time. This validates the 
procedure of energy minimization in relaxation theory. 

Magnetic helicity is equivalent to electron helicity H e 
for massless electrons: 



h! = H e = 



J A ■ V x Ad 3 r, 



(7) 



where A is vector potential, b = V x A. Magnetic he- 
licity is an ideal invariant; it changes only in presence 
of resistivity. For our study it is important that mag- 
netic helicity is more robust than the energy [1-3], i.e., 
its value is approximately constant in time and equal to 
its initial value, hi = K . 

As it has been mentioned above the ion helicity 



R = ar, t= —t, V = V A v, B = B Q b, P = pV\ V , 

VA 

where the Alfven velocity Va and the characteristic mag- 
netic field Bq are defined as 



Hi = / (A + <5 4 v) ■ V x (A + S^)d 3 r 



(8) 



is not a good invariant in Hall MHD. To illustrate this 
more precisely we introduce generalized cross helicity as 



(3) 



v (b + |v x v)d 3 r, 



(9) 
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which is related to ion hclicity by the equation: 



Hi — H„ 



26ihn 



i [ (v x A) • dS, 
Js 



where the last integral is taken over the cylindrical sur- 
face of the plasma. Though the time derivative of the 
integrand in Eq. (9) is reduced to full divergence, the 
generalized cross helicity (and, therefore, the ion helic- 
ity) is not an ideal invariant in Hall MHD, since 

where u = V x v is the fluid vorticity. In order for \i2 
to be conserved, an extra boundary condition must be 
imposed, which is u n = [9]. However, this boundary 
condition overspecifies the solution to Eqs. (1) and (2). 
Thus, the ruggedness of the ion helicity should not be 
assumed in Hall relaxation. 

If 5i = 0, Eq. (9) defines the MHD cross helicity: 



brf 3 i 



mi 



This quantity is an ideal invariant in incompressible 
MHD. Moreover, cross helicity (11) is also a rugged in- 
variant in MHD relaxation, which is confirmed by nu- 
merical simulations [8] . In our MHD analysis we assume 
that the cross helicity is constant and equal to its initial 
value, hi = M. 

The perfectly conducting cylindrical shell serves as a 
conserver of axial magnetic flux: 



h,= 



b z rdrdip. 



(12) 



o o 



This is the only true invariant in the system, which is 
conserved in both ideal and dissipative plasmas. Based 
on normalization (3), /13 = tt. 

The symmetry of the pinch configuration yields two 
more ideal invariants, the axial and angular momenta: 



hn = 



,d 3 l 



(13) 
(14) 



The initial values of these invariants can always be at- 
tributed to some uniform flow with velocity u in z- 
direction and a rigid rotation with angular velocity f2 
in (^-direction, i.e., /14 = nlu and = 7rZO/2, where 
I = L/a is the normalized length of the cylinder. 

We look for a minimum of energy (6) subject to con- 
straints given by global invariants (7), (11)-(14). The 
cross helicity (11) is used as a constraint only in stan- 
dard MHD but not in Hall MHD. The Lagrange multi- 
pliers method gives us a conditional extrcmum of E (it is 



minimum since E is positive definite, and, therefore, the 
resulting equilibrium is ideally stable): 



vo = -A*2t>o - /i 4 e 2 - n 5 re v , 

(1 - /uf)V x b + 2^ib = 2fi 2 [i 5 e z , 

hi = K, hi = AI, /13 = 7r, hi = irlu, /15 



(15) 
(16) 

(17) 



Note that the amplitude of the magnetic field b and 
the Lagrange multipliers Hi, (12, M4, ^5 are not arbitrary, 
they are determined by the constraints (17). Therefore, 
the relaxed state depends only on the initial values of the 
invariants. 

Eqs. (15)-(17) describe single-fluid MHD relaxed state. 
The Hall MHD relaxed state is obtained from these equa- 
tions by setting [ii = and ignoring cross helicity con- 
straint hi = M. The axisymmetric Hall state is 



An = 



1 



bn = 



v 



2Ji(A) 
A 

2Ji(A) 
Qre, n 4 



h(Ar) - J (A)e 3 
h(Ar), 



(18) 

(19) 
(20) 



where A = 
h(Ar) 



-2/ii, 



J (Xr)e z , V x h = Ah, (21) 



and Jo and J\ are Bessel functions of the first kind. For 
unique definition of vector potential An we take a gauge 
invariance condition [16] 



dz = 0. 



(22) 



Solution given by Eqs. (18)-(20) is in agreement with 
Ref. [13]; it corresponds to a force-free Taylor state [1-3] 
for magnetic field with addition of rigid rotation in az- 
imuthal direction or/and uniform flow in axial direction. 
Parameter A is determined from the magnetic helicity: 



K 



K 

Vi 



= 7 1- 



Jo(A)J 2 (A) 
J2(A) 



(23) 



Here and below we use tilde to denote the same quantity 
divided by irl, i.e., the normalized density of the quantity. 
For values of magnetic helicity larger than K > 4.1 the 
minimum energy state becomes helically symmetric [17]. 

Next we consider the single-fluid MHD relaxed state. 
In this case a solution to system (15)-(17) is 



A C 



h(Ar) - Ji(X)re v - J (A)e ; 



7^,(24) 



b = C I h(Ar) - ^pe : 



v = ~lxiC h(Ar) — re ¥ 



fire,. 



2Ji(A) 



(25) 
(26) 
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where h is given by Eq. (21), A and the amplitude C are 



A = 



2/ii 



C = 



[A(l - (4) - 2fyt 2 ] 
2[J 1 (A)+/i|J 3 (A)]' 



(27) 



Substituting these expressions into Eqs. (7) and (11) we 
obtain the magnetic helicity and the cross heicity: 



K 



2CJ 2 (X) 



A 



(28) 



2C 2 



J 2 (A)-2J (A)J 2 (A)- J 2 (A) 



M = u 



2ncj 2 {\) 

A 



(29) 



M2 C 2 2 J 2 (A) - 3 Jo (A) J 2 (A) - J 2 2 (A) 



8J 2 2 (A) 
A 2 



Eqs. (28) and (29) are used to find [i 2 and A through the 
initial values of helicities, K and M. In the limit /i 2 — >• 0, 
Eqs. (24)-(26) and (28) become Eqs. (18)-(20) and (23), 
respectively. We also note here that the uniform axial 
flow u can be removed from all equations. This corre- 
sponds to transformation to the reference frame moving 
with velocity u in z direction: 



vo 



M* = M 



Without loss of generality we assume that u = 0. 

Our main results are summarized in Figs. 1 and 2. A 
sample of MHD relaxed state (25), (26) with non-zero 
plasma flow is presented in Fig. 1. Fig. 2 shows the 
F — diagram of the relaxed MHD states with different 
values of cross helicity M and angular momentum f2. The 
reversal parameter F and pinch are defined here as 



F 




D- 



R=a 



B v | R =a 
Bn 



J(p \ r—l ■ 



(30) 
(31) 



As follows from Fig. 2, the presence of the initial flow 
(non zero cross helicity or total angular momentum) in 
cylindrical plasma pinch affects the relaxed magnetic field 
significantly. This is due to the coupling of the velocity 
and the magnetic field that occurs in Eq. (16) through 
the term on the right hand side. Such coupling does not 
take place in the systems without axial symmetry where 
the angular momentum (14) is not conserved, e.g., in a 
periodic rectangular box. In this relaxed mag- 

netic field corresponds to a force-free Taylor state and it 
has the same structure in both standard and Hall MHD, 
independent of the initial flows and determined by the 
value of magnetic helicity only [14, 15]. 

In summary, we have found the axisymmetric mini- 
mum energy (relaxed) states of a cylindrical plasma pinch 
within the framework of both standard and Hall MHD. 
Like all variational theories of plasma relaxation, the 
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FIG. 1: MHD relaxed state with magnetic helicity K — 2, 
cross helicity M — 0.1 and angular momentum Q, — 0.1. Solid 
lines are components of magnetic field (left vertical axis), 
dashed lines are components of velocity (right vertical axis). 
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FIG. 2: Reversal parameter F vs. pinch <d for the relaxed 
MHD states with different values of cross helicity M and an- 
gular momentum fi. 



present calculation is silent as to the details of the dy- 
namics that are responsible for the relaxation process. 
The only requirement is that they preserve the robust in- 
variants assumed during the variational procedure. Re- 
laxation with Hall MHD does not preserve the gener- 
alized cross helicity, leads to states with uniform axial 
flow and constant angular velocity, and reproduces the 
well-known Bcssel-function model for the magnetic field 
profiles. Relaxation with single-fluid MHD robustly pre- 
serves the cross helicity and leads to flow profiles with 
more radial structure, but it also modifies the form of the 
relaxed magnetic field. In both cases the solutions pre- 
dict redistribution (transport) of the momentum caused 
by plasma relaxation in a closed system. This conclusion 
is qualitatively consistent with flow dynamics observed 
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in RFP experiment. At this time there is insufficient 
experimental information to distinguish between MHD 
and Hall MHD relaxation. Further insight in this regard 
requires more detailed experimental measurements, and 
full nonlinear Hall MHD computations that take into ac- 
count the driven (open) character of RFP system. 

The authors wish to thank C. Hegna and C. Sovinec 
for useful discussions. This work is supported by the 
National Science Foundation and by the U.S. Department 
of Energy under Grant DE-FG02ER-54868. 



[1] J. B. Taylor, Phys. Rev. Lett. 33, 1139 (1974) 

[2] J. B. Taylor, Rev. Mod. Phys. 58, 741 (1986) 

[3] S. Ortolani and D. D. Scfmack, Magnetohydrodynamics 

of Plasma Relaxation, (World Scientific Press, Singapore) 

1993 

[4] A. Kuritsyn et al, Phys. Plasmas 16, 055903 (2009) 



[5] K. Chiyoda, J. Phys. Soc. Japan 55, 1139 (1986) 

[6] W. Schuurman and M. P. Weenink, J. Phys. Soc. Japan 

55, 4574 (1986) 
[7] J. Shiraishi, S. Ohsaki, and Z. Yoshida, J. Plasma Fusion 

Res. SERIES, 6, 169 (2004) 
[8] R. Horiuchi and T. Sato, Phys. Fluids 31, 1142 (1988) 
[9] L. Turner, IEEE Transactions on Plasma Science PS-14, 

849 (1986) 

[10] S. M. Mahajan and Z. Yoshida, Phys. Rev. Lett. 81, 4863 
(1998) 

[11] L. C. Steinhauer, H. Yamada, and A. Ishida, Phys. Plas- 
mas 8, 4053 (2001) 
[12] L. C. Steinhauer, Phys. Plasmas 9, 3767 (2002) 
[13] S. Ohsaki and Z. Yoshida, Phys. Plasmas 12, 064505 
(2005) 

[14] R. Numata, Z. Yoshida, and T. Hayashi, Comp. Phys. 
Comm. 164, 291 (2004) 

[15] R. Numata, Z. Yoshida, and T. Hayashi, J. Plasma Fu- 
sion Res. SERIES, 6, 130 (2004) 

[16] A. Reiman, Phys. Fluids 24, 956 (1981) 

[17] A. Reiman, Phys. Fluids 23, 230 (1980) 



